Heat conduction: hyperbolic self-similar shock-waves in solids 
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Analytic solutions for cylindrical thermal waves in solid medium is given based on the nonlinear 
hyperbolic system of heat flux relaxation and energy conservation equations. The Fourier-Cattaneo 
phenomenological law is generalized where the relaxation time and heat propagation coefficient 
have a general power law temperature dependence. From such laws one cannot form a second 
order parabolic or telegraph-type equation. We consider the original non-linear hyperbolic system 
itself with the self-similar Ansatz for the temperature distribution and for the heat flux. As results 
continuous and shock-wave solutions are presented. For physical establishment numerous materials 
with various temperature dependent heat conduction coefficients are mentioned. 
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In contemporary heat transport theory (ever since Maxwell's paper 1]) it is widely accepted in the literature that 
only for stationary and weakly non-stationary temperature fields the constitutive equation assumes that a temperature 
' ' gradient VT instantaneously produces heat flux q according to the Fourier law 

S 1 ' q(x,t) = -«VT(x,t). (1) 

Combining this equation with the energy conservation law the usual parabolic heat conduction equation is given. Heat 
i-^h , conduction mechanisms can be classified via the temperature dependence of the coefficient k ~ T u . There are three 
different cases of thermal conductivity, normal heat conduction which obeys the Fourier law (y = 0), slow (y > 0) 
and fast heat conduction — 2 < v < 0. 

In plasma physics if the temperature range is between 10 5 K and 10 8 K then the coefficient of the heat conductivity 
k depends on the temperature and density of the material. It is usually assumed to have a power dependence 
■ k = kqT^v 1 * where v = 1/p is the specific volume the coefficient «o and the exponents v, depend on the heat 
conduction mechanism 0. With radiation heat conduction one has 4 < < 6, 1 < fj, < 2; with electron heat 
conduction and fully ionized plasma v = 5/2, fj, — 0. For magnetically confined non-neutral plasma the classical heat 
' conduction coefficient is the following [| k k -j=ln[c2T 3 ^ 2 ]. Parabolic thermal wave theory is based on this approach 

0, Q ■ In plasmas heat conduction is strongly coupled to flow properties which we will not consider in the following. 

' The linear parabolic theory predicts infinite speed of propagation which is known as the "paradox of heat conduction" 
f*^S , (PHC). The following two theories resolve this contradiction. 

' However, if the time scale of local temperature variation is very small, Eq. ([T]) is replaced by 
i— I ' 

q(x,t + r) = -KVT(x,i) (2) 

• rH ' 

where r is called the thermal relaxation time. This is a thermodynamic property of the materials which was determined 
experimentally for large number of materials. Although r turns out to be very small in many instances e.g. is of 
order of picoseconds for most metals, there are several materials where this is not the case, most notably sand (21 s), 
H acid (25 s), NaHCOa (29 s), and biological tissue (1-100 s) @. 

Unlike the Fourier's heat conduction law, this constitutive equation is non-local in time. The desired local character 
can be restored with the Taylor expansion of q by time which is usually truncated at the first order namely 

q(x, i ) + r^M = - K VT(x,f). (3) 

This is the well-known Cattaneo heat conduction law Q the second term on the left hand side is known as the 
"thermal intertia". (Unfortunately, this form is not Galilean invariant, and gives a paradoxial results if the media is 
in motion, this problem was eliminated in by [7[.) Combining this constitutive equation with the energy conservation 
yields the hyperbolic telegraph heat conduction equation where r and k are constants. Hyperbolic equations usually 
ensure finite propagation velocity. Unfortunately, telegraph equations has no self-similar solutions which would be a 
desirable physical property. In the work of [§[ a non-autonomous telegraph-type heat conduction equation is presented 
with self-similar non-oscillatingcompactly supported solutions. A review with a large number of physical models of 
heat waves can be found in [5, 9]. A recent work on the speed of heat waves was published by [Io| . 
Our starting point is the following 
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qt = T r , (4) 

T T 

c T t = -q r . (5) 

r 

The first equation of the system is the generalized Fourier-Cattaneo heat conduction law and the second one is the 
energy conservation condition for the radial coordinate. The heat flux q = q(r,t) and the temperature dependence 
T = T(r, t) have radial coordinate and time dependence. The subscripts r and t notate the partial derivatives with 
respect to the radial coordinate and the time, respectively. (From now on we investigate the radial coordinate of a 
cylindrical symmetric problem as spatial dependence.) The parameter cq = pc where p is the mass density and c is 
the specific heat. Second order effects such as compressibility are neglected (p and are c constants during the process). 

In the following we shall suppose that the heat conduction coefficient and the thermal relaxation depend on tem- 
perature on the following way 

K — kqT u , t = t T-'. (6) 

The kq and tq are real numbers with the proper physical dimensions. Now our dimensionless system reads 



qt = -T e q-T e+ul T r , (7) 

T t = -<&•--• (8) 
r 

There are various phenomenological heat conduction laws available for all kind of solids, without completeness we 
mention some well-known examples. For pure metals according to |l2j (Page 275 Eq. 27.3) the Wiedemann- Franz 
low the thermal conductivity is proportional with the electrical conductivity a times the temperature k — aLT. 
The proportionality constant L is the so called Lorentz number with the approximate numerical value of 2.44 x 
10~ 8 WflK ~ 2 . For exact numerical data for various metals see [lj]. The relaxation time r is proportional to the heat 
conduction coefficient divided by the temperature. For metals with impurities the thermal resistivity (inverse of the 
thermal conductivity) is = AT 2 + BT^ 1 where A and B can be obtained from microscopic calculation based on 
quantum mechanics [12| (Page 297 Eq. 40.11). 

A hard-sphere model for dense fluids from [lj] derives a relation where the heat flux q(x,t) — aVT(x,t) + q 2 (x,t) 
which certainly means a non-linear heat propagation process. For the heat conduction in nanofluid suspensions 
[IH derives the k ss c/{T% — T{) law with additional time dependence. Another exotic and very promising new 
materials are the carbon nanotubes which have exotic heat conduction properties. Small et al. [lq performed heat 
conductivity measurements and found that at low temperatures there are two distinct regimes k(T) ~ T 25 (T < 
50K) and K(t) ~ T 2 (50 < T < 150isT). Beyond this regime there is a deviance from this quadratic temperature 
dependence and the maximum n value lies at 320 K. Above this value - at large temperatures - there is a k(T) ~ 1/T 
dependence according to [T3] - Additional nanoscale systems (like silicon films, or multiwall carbon nanotubes) have 
exotic temperature dependent heat conduction coefficients as well, for more see (l8j . For encased graphene the heat 
conduction coefficient is k ~ T@ where 1.5 < j3 < 2 at low temperature (T < 150K) [Hj]. A recent review of thermal 
properties of graphene and nanostructured carbon materials can be found in [20j . 

Our model is presented to describe the heat conduction of any kind of solid state without additional restrictions, 
therefore room or even higher temperature can be considered with large negative lj exponents. 

Even from these examples we can see that it has a need to investigate the general heat conduction problem, where 
the coefficients have general power law dependence. 

We look for the solutions of (|7I8[) in the most general self-similar form 

T = t- a f( V ) 7 q = t- s g(r 1 ). (9) 

For a better transparency in the following we introduce a new variable r) = where a, (3, S are all real numbers. 

The similarity exponents a, S and /3 are of primary physical importance since a, S represents the rate of decay of 
the magnitude T or q, while /3 is the rate of spread (or contraction if /3 < ) of the space distribution as time goes 
on. Self-similar solutions exclude the existence of any single time scale in the investigated system. 

We substitute © into J?]) and (0). It can be checked that 

a — — 7~p P= ot 1 | TT ' ^ oM"L e=ui + l. (10) 
CJ + 1 2(bJ + 1) 2(lj + 1) 
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Then we can obtain the shape functions f and g the following ordinary differential equation (ODE) system 

5g + Pruf = 9f u+1 + f 2u+1 f\ (11) 

(w)' = P(v 2 fY (12) 

where prime means derivation with respect to rj. 

The first lucky moment is that (1121) relates f and g in a simple way 

9 = M (13) 

if the a = 2/3 universality relation is fulfilled. 

Note, that we can immediately read how the self-similar solutions of the temperature distribution T and the heat 
flux q depend on uj 

r = i=rr/(^— , g = twi 9 — , (14) 

The parameter dependence of the complete heat conduction coefficient and relaxation time can be expressed via uj as 
well 



«=«ot^r — — , r = K t- L r +i — — . (15) 

Vt2(c+1) J \£ 2(^ + 1) J 

Recall that uj > — 1. These are already very informative and useful relations to investigate the global properties of 
the solutions, note that such kind of analysis are available for large number of complex mechanical and flow problems 

Substituting these relations back to Eq. (jlip after some algebra we arrive at the following non-linear first-order 
ODE 



(pf _ ^+1) = Pl [r+ i {2p + 1}]< (16) 



Put y = rj 2 and x = f . With this notation eq. (16) becomes linear for y{x) (this is the second lucky moment of 
investigation) : 



dy = y(x)-4(uj + l) 2 x 2 " +1 

dx x[(co + \)x u] + 1 - uj - 2]' 1 ' 

Plainly, / = is a solution to eq. (16). If y{x) the solution of eq. (17) is strictly monotonic then so is the inverse 
function f — x and no discontinuity. However if y(x) is not monotonic on some interval (x\,X2) and has a turning 
point at xoe(xi,X2) then the inverse (/ = x) has sense on [0,2/(2;)] only. One sets / = for y > y(xo) and the 
discontinuity a y(xo) is apparent. The analytical investigation of the linear equation (17) is, in general easier than of 
eq (16). In some cases (for some ujs) one can have more explicit or almost explicit solutions. 

There are two examples: 
The first case is for uj = 0, (a = 1, (3 = 1/2, 5 = 3/2, e = 1). 

This example was studied by [ll| in some details. The corresponding ODE (17) reads y' = (y — Ax)/x(x — 2) which 
has a solution y = 8 + [(x — 2)/x] 1 ^ 2 [ci — 8ln(^/x + s/x — 2)] where ci is a constant. 

It is clear that must be x > 2 and y{x) is monotonic for x > 2 until xq where y = 0. This means that x{y) exists 
and monotonic on some interval [0, yo]i ^(2/0) — 2; for y > yo we have x(y) = so the discontinuity. For a better 
understanding Figure la presents the graph of solution of Eq. (17) through the point (3,0.5). The inverse of this 
function for x > 2 is shown on Fig. lb (the nonzero part). The solid line is a solution through the /(0) = 10.8 point. 
Figure 2 presents the shock-wave propagation of the temperature distribution T(r, t) for uj = 0. 

The second case is for uj = -1/2, (a = 2, /3 = 1, S = 2, e = 1/2). 
Now Eq. (17) takes the form of & = 2(y - l)/[x{^/x - 3)]. It can be checked that y = c 2 x- 2 ^{x^ 2 - 3) 4/3 is 
a solution for any ci > 0. Take C2 = 1. The function y{x) is monotonic on (0, 9),y(9) — 0. Returning to original 
variables we have / = 9/[(rj 2 ) 3 ^ 4 + l] 2 (which is plainly less than 9!) According to eq. (14) temperature and heat flux 
distributions are 

( r 3/2 _|_ £3/2)2 ' 9 ( r 3/2 _j_ £3/2)2 ' ^ ' 
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y = r| 
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x = f 



f(T!) 6 
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a) 
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FIG. 1: The direction field of a) Eq. (17) for cj 
a) y(3) = 0.5 and for b) /(0) = 10.8. 



and b) Eq. (16) for to = The solid line presents numerical solutions for 



These solutions are not discontinuous. Analytical and numerical calculus suggest that oj = —1/2 is a critical exponent: 
for— 1<cj<— 1/2 the solutions are continuous, for cj > — 1/2 the shocks always appear. 
In summary 

We presented a hyperbolic model for heat conduction in solids where the relaxation time and heat conduction coefficient 
is a power law function of time. There are basically two different regimes available for different power laws. For - 
1 < uj < — 1/2 the solutions are continuous for all positive time and radial coordinate, for cu > —1/2 the solutions 
are only continuous on a finite and closed [0 : r/o] interval and have a finite jump at the the endpoint 770 ■ As physical 
interpretation numerous materials and solid state systems mentioned with temperature dependent heat conduction 
coefficients. 

The paper is dedicated to Annabella Barna who was born on 20th of December 2011. 
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FIG. 2: The shock-wave propagation of the temperature distribution of T(r,t) for ui — 
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